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The effect of a constant transverse magnetic field is analyzed for the propagation of VLF 
electromagnetic waves about the earth. These waves are considered to be propagating by 
modes in a parallel plate waveguide. The lower boundary of the guide is considered to be a 
perfect conductor, while the upper boundary is assumed to be a sharply defined, semi-infinite, 
homogeneous plasma with a constant magnetic field applied. The source of these waves is 
an infinite number of short, horizontal, cophased dipoles, uniformly spaced parallel to the 
constant magnetic field vector. Admittance relations are derived for the upper boundary of 
the guide by considering the waves to be of grazing incidence. It is found that these admit- 
tance parameters depend on the direction of propagation. These admittance parameters are 
plotted versus frequency for various heights. Expressions for eigenvalues and the eigen- 
functions are found as functions of the admittance. By employing a two dimensional 
Green's function, the amplitudes of the various modes due to the dipole source are evaluated. 
It is noted that the propagation constants differ depending on the direction of propagation, 
thus offering an explanation of the cast-west effect of VLF transmission. A numerical 
example is calculated and field strength versus distance values are found to correspond to 
some experimental results. 



1. Introduction 

The mode theory of VLF radio transmission has 
evolved from the study of electromagnetic wave 
propagation between the earth and the ionosphere. 
Wait [1] 2 and others [2, 3] have solved the VLF 
electromagnetic problem using mode theory and as- 
suming a sharply bounded ionosphere, while neglect- 
ing the earth's magnetic field. Friedman [4, 5] has 
made a modal analysis which differ in metliod from 
that of Wait, but arrives at similar conclusions. 
Friedman's analysis utilized an upper boundary of 
scalar impedance Z. 

All of these solutions have the same propagation 
characteristics regardless of the direction of j^ropaga- 
tion. However, Crombie [6] has noted evidence that 
VIjF propagation differs depending on whether the 
propagation is from east to west, or from west to 
east. He also indicated that this so-called '^east- 
west effect'' may be partly due to the earth's mag- 
netic field. Recently, Wait [7, 8] and Crombie 
[6, 9] have shown analytically that the rate of at- 
tenuation for east-west propagation behaves as had 
been indicated by the latter author. 

This paper attempts to show analytically the east- 
west effect on the basis of mode theory and iono- 
sphere-freespace boundary admittance. Because of 
the presence of the earth's magnetic field, the admit- 
tance of the sharply bounded ionosphere is not a 
scalar, but a tensor. Further, because of these mag- 
netic effects, the tensor depends upon the direction 
of propagation. 
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The solution of the radiation from a dipole in this 
parallel plate waveguide with anisotropic boundary 
conditions is, in general, difficult. However, one 
may solve an analogous two-dimensional problem. 
This simplifies the arithmetic while maintaining the 
qualitative structure of the three-dimensional case. 

The two-dimensional problem corresponds to the 
source of electromagnetic waves being an infinite 
number of horizontal, cophased, short dipoles, uni- 
formly spaced parallel to the direction to the earth's 
magnetic field. This situation is displayed in figure 1. 

In the sections that follow, the radiation from 
these sources are examined in detail. A numerical 
example is shown also. 
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Figure 1. Cophased sources for TM waves. 
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2. Dielectric Tensor of Magnetically In- 
fluenced Plasma 

The permittivity of an ionized gas may be deter- 
mined by examining the ion motion under the influ- 
ence of electromagnetic waves. Certain restrictions 
must be placed on this gas if the evaluation of the 
equation of motion is to remain simple. 

For the purpose of this analysis, the gas will be 
considered to be a homogeneous plasma. That is, 
the gas is uniform throughout and the number of 
positive charges is equal to the number of negative 
charges. The plasma exists in a state where colli- 
sions are relatively few. Since the positive and neu- 
tral particles have more mass than the electrons, 
their mobilities are correspondingly smaller and it 
shall be considered that their velocity is zero. Fur- 
ther, the number of electrons per cubic meter iV, 
does not vary appreciably with time. Thus the be- 
havior of a volume element of gas is the same as that 
of an ^'average^' electron obeying the laws of classical 
dynamics. These are the assumptions of the '^Mag- 
neto-Ionic Theory'' [10]. 

In the special case considered in figure 1, where a 
magnetic field is considered perpendicular to the di- 
rection of propagation, the ionosphere's relative di- 
electric tensor is given by, [11] 
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In expressions 2.2, 2.3 and 2.4 cop is the plasma fre- 
quenc}^, v is the collision frequency, cu^ is the cyclotron 
frequency, and a? is the frequency of the electro- 
magnetic wave propagating in the plasma region. 

Using 2.1, one is able to solve for the electromag- 
netic fields in the plasma region. The solution of the 
fields is required to examine the conditions at the 
plasma-freespace interface. 

3. Admittance Tensor of Plasma at 
Interface 

In order to solve for the electromagnetic fields for 
the problem described by section 1, one must first 
consider the restrictions placed upon the problem 
by the particular geometry. Since the sources are 
considered to be uniformly distributed along the y- 
direction axis, there is no ^/-variation for the fields. 



Further, the nature of these sources produces only 
transverse-magnetic waves. Considering these facts 
and eq 2.1, Maxwell's equations in the plasma be- 
come, 
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Solving eqs 3.1 in terms of the longitudinal electric 
field to obtain the following equation 
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The solution to this equation is of the form 

E,=Kie-'"e^'^' (3.3) 



where p and 7 are, in general, complex and the plus 
or the minus sign designates a wave traveling in the 
negative 2-direction or positive z-direction respec- 
tively. 

Substituting eq 3.3. into 3.2 one obtains the charac- 
teristic equation of 3.2 
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In the free-space region. Maxwell's equations may 
be solved for the transverse-magnetic case in terms 
of the longitudinal electric field yielding the following 
homogeneous wave equation 



The solution for the homogeneous wave equation 
depends upon the boundary conditions, and since the 
tangential electric field must be zero at a perfect 
conductor, the solution is of the form, 



Ez=K2 sin q^xe"^ 



(3.6) 



where q and y are complex. 

Equation 3.6 has a corresponding characteristic 
equation given by, 



-g'+y'+k'=0. 



(3.7) 
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From eqs ,S.l and 3.)^ the fields in tlie plasma region 
are related by the following relations 
jor the -f z-direction, 
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for the — z-direction 
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(cont'd) 



From the two tangential fields Ez and Tly^ an admit- 
tance relationship may be defined as 
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From the two expressions of eqs 3.8 and relation 3.9, 
one may obtain tw^o expressions for the admittance 
of the i^lasma boundary. 
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where Yg-j^ and Y^-e are the relative admittance 
parameters for east-west and west-east propagation, 
respectively. 

If grazing incidence is considered, upon the plasma- 
freespace interface, it is found that y-^jk, where the 
sign has already been determined by the form of eq 
3.6. This corresponds to the far-field situation. If 
this condition were not imposed, the admittance 
would be a function of the free-space eigenvalue g. 
Thus one obtains the following two expressions for 
the admittance of the plasma boundary, 
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The following restriction must be placed upon p, 

Re (i))>0, Im (p)>0. (3.12) 

With these restrictions in mind, the final expres- 
sions lor the admittances in terms of the constituents 
of equation 2.1 are 



Y =" 

•*• e-w 



j(4-v') 






(3.13) 



A numerical solution of these admittance parame- 
ters for height versus frequency is shown in figures 2 
and 3. The values of collision frequency and elec- 
tron density have been obtained from papers by M. 
Nicolet [12] and H. E. Houston [13], respectively. 
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Figure 2. Real part of admittance (Y12). 
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Figure 3. Imaginary part of admittance (Y12). 
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4. Eigenvalues and Propagation Constant 

In order to solve for a general expression for the 
eigenvalues in terms of the admittance matrix for all 
waves, one must also solve MaxwelPs equations for 
the transverse-electric waves. This solution, in terms 
of the longitudinal magnetic field Hz, is given as 



Hz=YoKs€'^''' cos qx. 



(4.1) 



Uecalling eq 3.6 and Maxwell's equations for TM 
and TE one can obtain, 
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Relating these four fields at x=a by the admittance 
matrix F, 
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Substituting the solutions 3.6, 4.1, 4.2, and 4.3 into 
4.4 and combining terms 
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For 4.5 to be true, the determinant 
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Evaluating the determinant 4.6, the following general 
expression for the eigenvalues is obtained in terms of 
admittance parameters. 
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where A^ = Fn F22 — F12 F21 . 

For the special parallel plate case under considera- 
tion, the only admittance parameter needed is Yu. 
Thus eq 4.7 becomes 
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Rewriting relationship 4.8, q=qn, n=l, 2, 3 . 
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If the term g^/fc F12 is large as is the case for the 
'^electric'' wall, the following approximation for the 
eigenvalue g„ may be obtained. 
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For the '^zeroth mode^' case, that is 71= 0, 
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and if n is greater than zero, qn may be approximated 
ut in the right-hand side of 4.10 yielding, 
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The propagation constant 
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may be obtained from the results of 4.11 and 4.12. 
Eestrictions on 4.13 are that 



Re(Yj>0,Im(7j>0; 



(4.14) 



when 4.13 is evaluated for n=0, the attenuation 
constant, Re(yn)=(Xn is very small, and for 7^>2, 
the phase constant, Im (t«)=/3«, is small represent- 
ing evanescent modes. 

One can see that if Fi2-^ oo , as in the case of a 
perfect conductor, expression 4.12 reduce to the 
familiar result, nw/a. Expressions 4.11 and 4.12 
are used, along with the data of Houston and Nicolet, 
to plot mode decibel loss versus distance on figure 4 
for n=Oj 1,2. The numerical values for the admit- 
tance used were 



Ye-w=22.5-37.5j, 
Yw-E=5Q.5-41.5j. 



(4.15) 
(4.16) 



This corresponds to a radian frequency of co=2X 10* 
and a height of 90 km. 

A synopsis of experimental results by Crombie [6] 
indicates that the ratio of electric field strengths 
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which is certainly true for the numerical case cal- 
culated and graphed on figure 4. 
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Figure 4. Attenuation versus distance. 



5. Mode Solution — Source Present 

With the source of current nioineiit IL eq 3.5 
becomes, 

(^2+£.+!c')KH'^f^JLd{x-b)5iz) (5.1) 

where d (^x—b) and 5 (i) are Dmic ^'deUa-f unctions. '^ 
This problem is analogous to the Green's function 
problem if one relates the longitudinal electric field 
to the Green's function as E2=^—jcjoixoILG, the cor- 
responding Green's function problem may be solved, 

{^2+^2+k')G=-S(x-b)S{z). (5.2) 

The boundary conditions for the Green's function 
are the same as for the electric field E, 



(G' + F(?)U,=0. 
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It has been shown in section 4 that the solution 
to the homogeneous wave equation for the longi- 
tudinal electric field and hence the Green's function 
is of the form 



6^=X;^A€-^"'-'^sin q,x. 
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Thus if 5.2 is true we must be able to expand the 
delta-function in terms of eigenvectors, [14, 15] 
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Multiplying through by sin q^n^x and integrating over 
the range of x, an expression for Bn 
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where r^ is often called the normalization factor. 
Integrating 5.2 over the source with respect to z 
yields 

Vn sin qnb 

and 



^n " 



(5.8) 



G=J2" 



2?„ 



r' 7«(2(Z„a— sill 2q„a) 



e ''"l^l sin 2„6 sin g^x. 



(5.9) 



Therefore the expression for the electric field is given 

by 
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The other fields may be evaluated for Maxwell's 
ecjuations for the TM case. 

The expression for the vertical component of the 
electric field with a ^^normahzed^' current is given as 



E.= ^n 



sin q^b 



T' 2g^a— sin 2g„a 



g-7„|zi (.()g g^^^ (5.11) 



If the transmitter is in an airplane at a height of 
a/4 and the receiver is on the ground, 5.11 becomes, 
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For tlie first three modes 5.12 is of the form 
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The magnitude of the field strength at the receiver 
is given by the expression, 
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where \Cn\ and Cn are the magnitude and argument 
of Cn^ respectively. 

Equation 5.15 is plotted on figure 5. These 
results depend not only on attenuation, but also on 
field strength magnitude which in turn depends on 
the direction of propagation and the height of the 
transmitting antenna. This corresponds qualita- 
tively with the far field experimental results men- 
tioned by Crombie [6]. 
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Figure 5. Field strength versus distance. 



6. Conclusion 

The effect on east-west, west-east propagation of 
electromagnetic transverse magnetic waves across 
the earth^s magnetic field has been analyzed as a 
boundary value problem. Certain idealizations were 
made to solve this problem. Admittance relation- 
ships were derived from field relations existing in a 
homogeneous anisotropic plasma, and it was found 
that the admittance parameters depended upon the 
direction of propagation. These parameters were 
plotted versus frequency for various heights. Ex- 
pressions for the eigenvalues were derived. The 
attenuation over distance was plotted for a special 
case. It was found that the ratio of the electric 
field for the individual modes of west-east/east-west 
propagation was greater than one. This corresponds 



to some experimental data. A two dimensional 
Green's function was derived and a numerical calcu- 
lation of the corresponding electric field strength was 
plotted versus distance. 
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